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o ■ 

We revisit the problem of the spin-Peierls instability in a one dimensional spin-^ chain coupled 
£SJ ' to phonons. The phonons are treated within the mean field approximation. We use bosonization 

! i _ techniques to describe the gapped spin chain and then use the Thermodynamic Bethe Ansatz to 

O ■ obtain quantitative results for the thermodynamics of the spin-Peierls system in a whole range of 

temperature. This allows us to predict the behavior of the specific heat and the magnetic suscepti- 
bility in the entire dimerized phase. We study the effect of small magnetic fields on the transition. 
Moreover, we obtain the parameters of the Landau- Ginzburg theory describing this continuous phase 
transition near the critical point. 

PACS numbers: 75.10.Pq 75.40.Cx 75.50.Ee 

; 

The spin-Peierls instability^ is the magnetic analogue of the Peierls instability of a one-dimensional metal^ In 
the spin-Peierls instability, an antiferromagnetic spin-1/2 chain coupled to optical phonons develops a spin gap via a 
static deformation (or dimerization) of the lattice at zero temperature. In this dimerized phase, the gain in magnetic 
, energy resulting from the formation of the spin gap outweighs the loss of elastic energy due to the static deformation. 
For a system consisting of an array of spin chains coupled to two or three dimensional phonons, the dimerized phase 
can persist for temperatures < T < T$p, where T$p is the spin-Peierls transition temperature. For T > T$p, 
■ . the chains are undistorted and and one recovers gapless spin excitations. The phase transition at T = T$p between 
the dimerized and the uniform state is a second-order phase transition and has been observed in a host of quasi-one 
dimensional organic materials such as TTFCuS 4 C4(CF 3 )4 and TTFAuS 4 C4(CF 3 )4 (also known as TTFCuBDT and 
TTFAuBDT)4, MEM-(TCNQ)A(TMTTF) 2 PF 6 , (TMTTF) 2 AsFg&L£, and (BCPTTF) 2 PF 6 £ The discovery of the 
l—— inorganic material^ CuGeOa, spurred further activity in this domain as this system is more convenient for neutron 

■ scattering studies. 

From the theory point of view, most of the treatments consider the phonons as statici^iAi This mean field treatment 

■ of the phonons is expected to work when the phonon frequency can be neglected (adiabatic limit) compared to the 
\Q ' spin gap. Such an approach is thus better suited to the softer organic materials than to the inorganic compound 

: CuGeOa where the frequency of the phonon driving the transition is not small compared to the spin gap^ However, a 
common feature of the spin-Peierls transition in all the spin-Peierls compounds, is that some data indicate a BCS type 
mean field behavior of the thermodynamic quantities near the transition^ For instance, a BCS type relationship^, 
A/ksTsp = 1.76, between the zero temperature spin gap and the spin Peierls transition temperature Tsp has been 
observed in CuGeOa, and it was used to argue that the transition in this material could also be described within 
mean field theory. However, the exact nature of the transition in CuGe03 is still disputed^ In particular, no 
phonon softening was observed near the transition^ in disagreement with the mean-field scenarioiSiii The absence 

1 i of phonon softening could be attributed to the high frequency of the phonons coupled to the spin excitations^ 

Other discrepancies with the mean-field scenario are discussed in Ref. [RJ Nevertheless, despite these deviations, 
phcnomenological Landau-Ginzburg theories can be used to some extent to fit the critical behavior of CuGe03. 15 
Therefore, it is important to develop a more quantitative description of the mean field theory of the spin-Peierls 
transition, particularly in the gapped phase in order to have a more reliable comparison of the predictions of the 
^ ■ mean-field scenario with experimental data on the spin Peierls materials. 

k^l ' In the first theoretical approaches to the spin-Peierls transitions^, the spin chain was mapped onto a model of 
interacting one-dimensional spinless fermions by the Jordan- Wigner transformation^ and the interactions between the 
fermions were either neglected 1 or treated in the Hartree-Fock approximation^ Later, in Ref. [TJ, though the phonons 
were still treated at the mean field level, the spin chain was described using bosonization which correctly describes the 
quantum critical behavior of the pure spin g chaini^SiSSiSii 2 ^ A linear response treatment of the spin-phonon coupling 
resulted in a much improved estimation of dependence of the transition temperature on the spin phonon interaction. 
Furthermore, a Landau Ginzburgexpansion was developed to study the vicinity of the transition. 

However, in contrast to Refs. [iJQ no prediction could be made for the thermodynamics in the dimerized phase. 
There are two reasons for this. First, in Ref. I 111 the dimerized phase is not described by a model of noninteracting 
fermions with a gap as in Ref. Q], but by a more complicated massive sine-Gordon models Second, in bosonization, 
although the expressions of the lattice operators in terms of sine-Gordon fields are known, the amplitudes in this 
expressions were unknown and have been determined quantitatively only recently^*2£*2& Since the thermodynamics 
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of the massive sine-Gordon model is now understood2L22i2i!iSL2ii2£, and the exact expression of the gap in the sine- 
Gordon theory is known™,, it is now possible to study the thermodynamics of the dimerized spin-Peierls phase within 
the mean field approximation, as well as study the zero temperature properties. 

In this paper, we will use the above developments to revisit the problem of the spin-Peierls transition in the adiabatic 
approximation for the phonons. Our methodology and results are also applicable to the chain mean field theory of 
quasi-one dimensional antiferromagnets since the latter presents a formal analogy to the theory of the spin-Peierls 
statei2^i25i2£i2L22i2£ The paper is organized as follows: In Sec. |U we present the model and its bosonized version. We 
obtain analytical results results for the spin-Peierls temperature T$p and the total energy and gap at zero temperature. 
In Sec.[n] the Thermodynamic Bethe Ansatz is used to study the thermodynamic properties of the dimerized chain at 
finite temperature. We obtain various results for the gap, the dimerization, the specific heat and the static magnetic 
susceptibility for an entire range of temperatures smaller than the bandwidth. We use these results to derive the 
Ginzburg Landau functional describing the mean field transition and study the behavior of the correlation length near 
the transition. 



I. MODEL 



Within a mean field treatment of the phonons, the full Hamiltonian describing the coupling of the lattice to the 
spin chain is given by 



§« + J(l + (-) n A(«))S n -S n+1 



(1) 



where (u) is the mean- field displacement along the chain and the Si are spin-i operators. K is the elastic constant for 
lattice deformations and A is a parameter related to the amplitude of the spin phonon coupling. For the sake of clarity, 
we introduce the dimensionless variable S — X(u), and the reduced elastic constant K = K/X 2 . To obtain the physics 
of the spin-Peierls transition, we need to evaluate the mean field displacement (it), or equivalently, the parameter 5 
in a self- consistent manner i.e., the value of S which minimizes the free energy at any given temperature T. In the 
spin-Peierls phase T < Tsp, S(T) ^ and S(T) = for all temperatures T > Tsp- A self-consistent evaluation of 
S(T) then permits a systematic calculation of various properties of the spin chain in a rather straightforward manner. 
Since we are mainly interested in the low energy/long wavelength physics of the spin Peierls system, we study the 
Hamiltonian Q in the continuum limit. In this limit, the continuum Hamiltonian reads: 



// = / dx^S 2 



where, k — K/a, a is the lattice spacing and H s is the continuum approximation of the spin Hamiltonian. 
standard bosonization techniquesi2*2ii22iaS, the continuum spin Hamiltonian is found to be : 



dx 
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Using 
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where, the fields </> and n are canonically conjugate to each other (i.e. [cf)(x),Il(x')] = i8(x — x')). The velocity 
of the bosonic excitations defined by the field cj> is u = \ Ja and g is an amplitude proportional to the exchange 
interaction J. The Hamiltonian J3J| is the well known sine-Gordon model with 1 = 2ir. In Eq. ©, we have omitted a 
marginally irrelevant term^i We will come back later to the effect of this term on the properties of the system. In the 
Hamiltonian a non-zero dimerization 5 induces the relevant operator, cos\/20 of dimension i. This results in a 
gap A ~ <5 2 / 3 , and a diminution of the ground state energyH^i E s (5) — E s (0) ~ — S 4 ^ 3 . For small 8, this reduction of 
magnetic energy compensates the loss of elastic energy in resulting in a dimerized state at T — 0. Until recently, 
the proportionality constant between g and J was unknown, thus preventing a quantitative estimation of the magnetic 
energy E(S) and hence the correct value of the spin gap A. Consequently, only exponents could be predicted from the 
above mean field description, and no prediction could be made for the thermodynamics of the system below the spin 
Peierls transition temperature. However, recent developments in integrable systems and bosonization, now permit a 
precise determination of the amplitudes in the continuum theory^25^ A correct mapping of the lattice spin model 
onto its bosonized version fixes the amplitude g in J3J : 
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Although the present paper focuses on the spin-Peierls system, we reiterate that the approach outlined below, is also 
applicable to the chain mean-field theory of quasi-one dimensional antiferromagnets 3 ^ 3 ^* 3 ^^ 3 ^ 3 ^ as long as marginal 
operators are neglected. In the case of the antiferromagnet, the magnetization m and the inverse of the interchain 
exchange term J J 1 play the role of the dimerization and the elastic constant respectively. We now analyze the full 
Hamiltonian J2J in certain limits. 



A. Zero temperature limit 

At zero temperature, the dimerization 5 is non-zero, resulting in a gap for spin excitations. As mentioned earlier, the 
precise mapping of the spin lattice model onto its continuum version, the sine-Gordon model, yields exact expressions 
for the gap and the total energy of the spin system^ In this model, the lowest energy excitation is a soliton 2 ^ and 
using 10}, its mass is given by 



M = AEilM 

0FT(2/3) 



r(3/4) 3 /Try/4 
r(l/4)7r 2 \2J ' 



2/3 

(5) 



Besides the soliton and the corresponding antisoliton excitations, the sine-Gordon model at (3g G = 2ir possess two 
other excitations, a light breather with a mass M and a heavy breather with a mass My/3. The soliton, the antisoliton 
and the light breather together form a SU(2) spin triplet while the heavy breather forms a SU(2) singleti^ 2 * 4 ^^ The 
gap to the lowest energy excitation or equivalently, the singlet-triplet gap, is 



A = —M ~ 1.723J5 2 / 3 (6) 
a 

A comparison of this predicted value with the real gap of the spin lattice system calculated numerically using the 
density matrix renormalization group was done in Ref. |4jj and a reasonably good accord was found. The knowledge 
of the soliton mass (JSJ also yields the ground state energy per spin of the dimerized spin chain 3 ^: 

E s (5) = t an \ - -0.2728 J<5 4/3 , (7) 

which is in reasonable agreement with numerics4^ To obtain the effective dimerization S at zero temperature, we need 
to minimize the total ground state energy per unit spin of the spin lattice system, E = K/25 2 + E S (S) which leads to 
the following results: 
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and we indeed have the 
However, the prefactors 



We note that the ratio J 2 j K can be identified with the coupling constant Xqp used in Re f . ITll 
same exponents for the dependence of the gap on the coupling constant Xcf as in Ref. I 111 ' 
in Eqs. © could not be obtained in Ref. [Ill 

The continuum approximation underlying our mean field theory is valid when the zero temperature correlation 
length is much larger than the lattice spacing i.e. when Ac J. Clearly, this requires that J -C K i.e., a sufficiently 
rigid lattice. This criterion leads to S(T = 0) <C 1. We note that relations similar to the ones in © have been 
previously obtained in the context of the chain mean field theory of quasi-lD antiferromagnetsrfS In reality, the 
results of ijHJl are slightly modified by the presence of a marginally irrelevant term in the continuum Hamiltonian 
for the spin system4i^£*^£ When the marginal interaction is taken into account, it is found that in the spin-Peierls 
case A ~ 5 2 / 3 | ln<5| -1 / 2 , whereas in the case of the antiferromagnet A ~ S 2 ^ 3 \ Ind] 1 ^ 6 i.e. the marginally irrelevant 
term frustrates the dimerization and favors antiferromagnetic ordering. The marginal interaction is eliminated in 
the J1-J2 chain at its critical point^i 4 ^ J 2 /Ji ~ 0.24. With an additional dimerization of the nearest neighbor 
exchange^ in this critical chain, the gap A = 1.76<5 2 / 3 . Moreover, in the absence of dimerization^, the spin velocity 
at the critical point is found to be 1.1936Jia. Generalizing the results obtained above, we find that the following 
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amplitude g = 0.8067r 2 Jia should be used in the bosonized Hamiltonian © in order to describe the J\ — J2 chain 
with J2/J1 = 0.2411. The resulting energy gain from dimerization is then E s /Ji = — 0.3745<5 4 / 3 and one obtains the 
following zero temperature results: 
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A = 0.879^- 
K 

(U.V, ( h 
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/3 

E = -0.093^ (9) 

Comparing and (jHJ , we see that the introduction of J2 results in a st rong enhancement of the zero temperature gap 
and of the dimerization, in agreement with a scenario proposed in Ref. 15 ll for the spin-Peierls transition in CuGe03. 

B. Transition temperature 

In this section, we redo the calculation of Ref. [ll| yielding the spin-Peierls transition temperature. For any tem- 
perature, a self-consistent treatment of the problem requires a calculation of the free energy as a function of the 
dimerization S taken as a variational parameter, followed by a minimization with respect to S. Unlike the zero 
temperature case, calculating the free energy for arbitrary temperatures requires the use of Thermodynamic Bethe 
Ansatz techniquesSLSS* 2 ^, and no closed analytic expressions can be obtained. However, to calculate the spin Peierls 
transition temperature, this full treatment is not required^! Indeed, close to the spin-Peierls transition, the order 
parameter 6 becomes small and a second order perturbation theory in 5, is sufficient to evaluate the leading behavior 
of the variational free energy. A straightforward perturbative development in S of the Matsubara imaginary time path 
integral gives the following expression for the free energy of the sine-Gordon model: 

r dx r iT £ + (10) 

611 ifiu V(2iro)V J-o* Jo Jcsh^-cos^ 

Note that the first term in this expression is just the free energy of a non-interacting Bose gas in one dimension. Using 
Eq. (8.12.4) in Ref. |52jto integrate over the space variable x, we obtain 

'•<™--s^-iG&)YM-""x)*- 

where the function P_i/2 is a Legendre function. A final integration over r using Eq. (8.14.16) in Ref. |S3 leads to: 
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The full mean field variational free energy is F MF (T,5) = \8 2 + F S (T,8) = CS 2 + o(5 2 ). When C > 0, which is 
obviously the case for high temperature, the mean field free energy has a minimum for 5 = and for C < 0, the 
energy is minimized by a state with non-zero dimerization. Therefore, the spin-Pcicrls transition temperature is fixed 
by the condition C — and using (|12|) . we obtain 

9 /7T\1/2J 2 „_„J : 



(I) ^ = 0.25342^ (13, 



27r 2 P(3/4) 4 

Note that the validity of the continuum description requires that T$p <C J. In Ref. ITTl the same dependence of 
Tsp on J 2 /K (up to the prefactor) was derived using an equivalent response function formalism. Comparing the 
two expressions, we observe that in Refiii, the transition temperature T$p — 1.01J 2 /K obtained there is a gross 
overestimate of T$p highlighting the importance of having correct amplitudes in the bosonized theory. Comparing 
Eqs. © and i|13|) . we note that the ratio: 

^2 = 2-47 (14, 

J-SP 
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is independent of the various coupling constants present in the theory. This ratio is in accord with values obtained 
by numerical studies of the spin-Peierls problem^ The existence of such an universal ratio is reminiscent of the BCS 
mean field theory for superconductivity^ where the ratio of the superconducting gap and transition temperature is 
approximately 1.76 . In fact, one can use the Jordan- Wigner transformation^ to map the Heisenberg spin chain onto 
a chain of interacting spinless fcrmions. Neglecting the interactions^, the resulting theory presents a formal similarity 
with the BCS theory^, which leads one to anticipate an universal ratio. We note, however, that the fact that the 
spinless fermions theory is strongly interac ting renormalizes the BCS ratio away from the non- interacting value 1.76. 
In particular, as already discussed in Ref. 1 131 the observation of a ratio of 1.76 between the zero temperature gap 
and the transition temperature in CuGeOs cannot be taken as an indication of adiabatic behavior in this compound. 
As discussed earlier, the results of (jHJ were obtained neglecting the logarithmic corrections induced by a marginally 
irrelevant interaction. These marginal interactions affect the dependence of the gap and the ground state energy 
on the dimerization, particularly for 8 <C 1 and at finite temperatures induce logarithmic corrections in response 
functions^ This inhibits a precise estimation of the BCS like ratio especially in systems with a small dimerization 
at low temperatures. 

For the next nearest neighbor chain with a critical coupling J2 C = 0.2411, where logarithmic corrections vanish, 
A = 1.5386Jf/K and T S p = 0.623J?/K. Note that these values respect the BCS relation JTI). In the light of the 
preceding discussion, it is interesting to note that a small change in the velocity and the coefficient of the sine-Gordon 
term, leads to a big change in the gap and the spin-Peierls temperature. This implies that the frustration engendered 
by J<i enhances fluctuations towards spontaneous dimerization, hence favoring the formation of the spin-Peierls state. 



C. Effect of logarithmic corrections 



We now discuss the effect of the marginally irrelevant operator cosy 80, neglected in the preceding sections. 
This operator is known to induce logarithmic corrections in the dimerization garj 4 ^ 4 ^ 4 ^ as well as in response 
functions. 48 i 55 i 56 i 57 This results in a modification of Eqs. © and (|13fl and consequently deviations from the BCS 
like ratio 114|) . Including these corrections, the gap at T — and the ground state energy are: 

A ~ J WsJJ^ (15) 

8 4 / 3 K 

Minimizing the ground state energy with respect to 8, one finds: 

«V»|M~-jp (17) 

For the transition temperature, it can be shown following Ref. [S3 that the susceptibility associated with the dimer- 
ization operator is corrected by a logarithmic factor so that: 



-3/2 
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With this result, the equation l|13|) is modified into: 



J \ 3/2 J* 



Ts P \ln—) ^ (19) 

We see from Eqs. ()17(l and l|19(l that the effect of logarithmic corrections is to decrease the spin-Peierls transition 
temperature and the zero temperature gap. In contrast, in the case of the Neel state, these logarithmic corrections 
enhance the transition temperature and the order parameteri^iiS^LiSiiS 

The equation l(T7jl leads to: 

/ J. \ 3/2 

(20) 
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resulting in a gap: 



^~^>—^i (21) 
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Solving for Tsp in Eq. I|19|l . and comparing with Eq. (|21[) . we find that to lowest order a BCS type relation holds. 
This relation however is obtained by retaining only the lowest order logarithmic corrections. 

II. THERMODYNAMIC BETHE ANSATZ MEAN FIELD THEORY 
A. Mean Field equations 

We now focus on the mean field theory for the spin lattice system at arbitrary temperatures. We use the Thermody- 
namic Bethe Ansatz (TB A) as a tool to evaluate the finite temperature free energy of the sine-Gordon Hamiltonian 10 • 
The TBA treatment of t he generic sine-Gordon model with a relevant term cos /3<fi has been formulated using the 
string hypothesis in Refs. 27 28 29. In general, this leads to an infinite number of coupled integral equations for the 
various pseudoenergies. However, at the so called reflectionless points^, defined by f3 2 = ^L, where n is an integer, 
the number of independent integral equations becomes finite^ Numerical methods can then be used to solve these 
integral equations and deduce the thermodynamics. The case of the dimerized spin chain with 1 — 2ir falls into this 
category and the TBA equations of Refs. B7l28ll2^ can be used to calculate the free energy. For generic values of (3 
away from the reflectionless points, the general formalism developed by Destri and de Vega^^i^ is more appropriate 
than the string approach. This latter method has been used successfully to obtain the thermodynamic properties of 
copper benzoate^ in a magnetic field. This approach can be used to study the thermodynamics of the generalized 
spin-Peierls transition 6! 2*£i or the antifcrromagnetic transition^ in systems of coupled spin ladders in a magnetic field. 
Before we embark on an application of the TBA method to the spin Peierls system, we note that in Refs. I27l28l29l 
the free energy is taken to be zero at zero temperature. However, since our reference state is the undimerized chain, 
we must add the zero temperature di merizatio n energy to the free energy calculated using the TBA of Refs. I2?ll28ll29l 
Using the approach outlined in Refs. I27l28l29l we find that in our case, the sine-Gordon free energy reads: 

FJT,5) = —— d^Acosh^ 31n(l + e- £1 W/ T ) + V31n(l + e- £2 W/ T ) " tan ^i!l_ ( 22 ) 
2nu J_ 00 L J a z 6 4 

where the pseudoenergies ei{9) and e 2 (9) are self-consistently determined by the following integral equations: 

ei(0) = Acosh<9 + — / dd'K 11 (e-9')ln{l + e~^ (e ">/ T ) + — / d9'K y2 {9 - 6 1 ) ln(l + e - £2 ^/ T ), 
2tt J^oo 2tt 

OTi r CO rp p CO 

e 2 (0) = AV3cosh0+^- / d9'K 12 (9-8')ln(l + e- ei( - e ' )/T ) + — / d8'K 22 (8 - 9') ln(l + e^ e ">/ T ), (23) 

2 ?r J-oo 2lr J-oo 



The integral kernels are given by: 



K 12 {9) 



2 sin f cosh 9 
sinh 2 9 + sin 2 | 
2 sin 5 cosh 9 2 cosh 9 sin ^ 

6 _|_ 2 

sinh 2 9 + sin 2 | sin 2 £ + sinh 2 9 



2 

K 22 {9) = 3K n (9) (24) 

The pseudoenergies ei,2 have a transparent physical interpretation^: t\{9) represents the dressed energy of the 
solitons and of the light breather (which have identical masses at the 1 = 2ir point) , whereas the pseudoenergy e 2 (9) 
represents the dressed energy of the heavy breather. In fact, because scattering is diagonal, Eqs. (|22|) - (|24|l can be 
easily re-derived using the approach outlined in Ref. I63L It is useful to recast the dimensionless free energy / — aF/J 
in terms of the scaled energies ti = £i/J and the reduced temperature T — T j J. 
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/ = — S 2 - — f dOMcoshe [3In(l + e" 6 "" 1 ^ ) + V3h(l + e - e »W/ f )l _ ^-M 2 (25) 
2J 2vr7 L y J 8\/3 

OTi /• CO rp pOC 

eUO) = -Mcosh6> + — / d6»' J ft: 11 (6l-6' , )ln(l + e~ el(e ' )/f, ) + — / dfl'iC^ - 6') ln(l + e -^')/?) ( 2 6) 

2 27T 27T J_ x 

O .'Ft /» CO rp pOO 

e 2 (6) = -Mv / 3cosh6»+ — / d6'K 12 (6 - 6') ln(l + e - ll{B ' )/f ) + — / d9'K 22 {6 - 6') ln(l + e ~ l2(e ' )/T \27) 

2 27T J_ 00 27T J_ 00 

To obtain the free energy for arbitrary temperatures we need to solve for a given soliton mass M, the pair of self 
consistent equations \'l'2l for the dispersions t\ and e 2 - Here, we use numerical methods 2 ^ to obtain the solutions for 
any arbitrary temperature. We use a simple iterative procedure to solve (|26JI - (j27(l numerically for various values of 
the dimerization S at a fixed temperature T. This provides us with the variational free energy energy for an entire 
range of 6 at fixed temperature. We then identify the value of 5 for which the free energy is a minimum. Such a 
procedure is repeated for various values of the temperature T thereby permitting us to obtain S(T). In particular, 
we find that the numerical estimate of Tgp coincides perfectly with the prediction of perturbation theory Eq. . 
This provides a first check of the validity of our TBA mean field solution. In Fig. [21 we plot our results for the mean 
field dimerization S(T) as a function of the reduced temperature. These results for S(T) are then used to obtain 
thermodynamic quantities. Figs. ^ an d ID show a plot the gap and specific heat as as function of the temperature. 
In the vicinity of Tgp, the gap vanishes as A a (Tsp — T) 1 / 3 and the specific heat jumps at T$p as expected for a 
second order transition. For low temperatures T <gC Tsp the specific heat is exponentially suppressed by the spin gap 
and for T >T$p, the specific heat is simply that of the pure Heisenberg chain: C v = ^-T — |j which is the same as 
that for a gas of free bosons. 



B. Law of corresponding states 

The mean-field theory leads to a law of corresponding states^ or equivalently, scaling forms for the free energy and 
associated quantities. We use the T = result S ~ (J/Kf/ 2 (cf. Eq. ©) to rewrite the finite temperature dimerization 
8(T) = {J/Kf/ 2 5{T). Inserting this in Eq. ©, and in Eqs. (H2J- and using H13JI . it is straightforward to see 
that the pseudocncrgies satisfy the scaling form ei(9) = Tspe(T/Tsp,0)- Consequently, this implies that the total 
free energy, gap and dimerization can be re-expressed as: 

F(T) = (-Pi (28) 



S(T) = l-^J 51 — ) (29) 



J J \T S p 

A(T) = T SP A (JL} (30) 

where the functions A, 5, f are universal functions of the scaled temperature. From Figs. ^ an d 121 we see that the 
numerical solutions for A, 6 do obey the above scaling form. From the expression for the free energy ()28|l . one easily 
obtains the following result for the specific heat: 

Tsp T ( T 



c - ^ -TttJ {%;) (31) 

As expected in a mean-field theory, there is a jump in the specific heat at the transition whose magnitude is given by 

The numerical solution of the mean field equation yields jsp = 1-39. Our result for the specific heat is shown in 
Fig. |3 A universal ratio of the specific heat jump to the specific heat above the critical temperature^ exists in 
the BCS theory, where AC v /C v (T+ p ) = 1.43. Here again, the value of the ratio AC v /C v (Tg P ) = 2.1 in the spin 
Peierls problem is different from the BCS ratio of 1.43 due to the strong interactions between the Jordan- Wigner 
pseudofermions. We note that in experiments on CuGeC>3, this ratio was found to be 1.5 or 1.6, which is close to 
the BCS prediction^ However, as in the case of the ratio of the gap to the transition temperature, this apparent 
agreement with the mean field description proves to be spurious as we have seen that the ratio should be near 2.1 in 
a spin isotropic material. 
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FIG. 1: The dimensionless scaling function A describing the law of corresponding states followed by the spin-Peierls gap. The 
universal ratio 2.47 is reached for T < 0AT SP . For T -> T S p the scaling function A ~ (1 — T/T S p) 1/3 . 
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FIG. 2: The dimensionless scaling function S describing the law of corresponding states followed by the spin-Peierls dimerization. 
The zero temperature value is reached for T < 0.4Tsp. For T — > Tsp the scaling function S ~ (1 — T/Tsp) 1 ^ 2 . 
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FIG. 3: The specific heat of the spin-Peierls problem in the mean field approximation. 
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C. Landau-Ginzburg expansion 



In this section, we use the results of the preced ing sections to obtain a simple Ginzburg Landau functional describing 
the vicinity of the spin-Peierls transition. In Ref.[Tl], it was shown that a soft Ising or 4 theory was enough to describe 
the vicinity of the transition. However, the coefficients, in particular, that of the quartic term could not be entirely 
calculated. Our formalism permits us to obtain the leading terms in the functional with the correct prefactors. This 
will be useful for more sophisticated treatments of the transition taking into account the fluctuations of the latticed 
or the role of solitons in the thermodynamics^ A Landau Ginzburg expansion^ of the variational free energy per 
unit length ll'L'l) in the vicinity of T$p gives 

F(T,<5) = |(T-T SP )«S 2 + |<5 4 (33) 

The law of corresponding states (|28|l leads to some constraints on the form of the expansion. Minimizing with respect 
to 5 yields: 

S 2 (T) = ^(T SP -T) 

2 

F(T) = -^-(T-T SP ) 2 (34) 

The law of corresponding states (|28(1 implies that p/q ~ T| P /J 3 and that p 2 /q ~ 1/J. Thus, we have p = c\J 2 /Tg P 
and q — C2J 5 /Tg P , where ci, Ci are dimensionless numbers. These predictions are in agreement with the ones obtained 
from the RG treatment in Ref. |6^. The dependence of p can also be verified by the perturbation theory of Sec. II Bl 
Although the precise value of q was not obtained in Ref. Illl it was shown using perturbation theory that the Landau- 
Ginzburg free ener gy h ad an expansion in powers of {J/T SP ) 1/2 A /T SP , where A = JS, Reporting this expansion 
in Eq. (5.5) of Ref. Illl leads precisely to the dependence of q on J and T$ P . Thus, the perturbative expansion of the 
free energy is fully consistent with the law of corresponding states. In terms of the dimensionless constants, 

ci / T SP \ 3 / T \ 



= _± (T-T SP ) 2 

4C2 J 

This also implies from Eq. (J5J that the spin-Peierls gap vanishes as ~ (1 — r/Tgp) 1 / 3 near the transition. This 
behavior is entirely consistent with our numerical results Eq. 134(1 . It now suffices to calculate the constants C\ and 
C2- A comparison with Eq. Q12JI. fixes c\ = 0.2534 and the value of C2 = 0.02276 is obtained by fitting the TBA mean 
field theory results for S(T) to $Q in the range 0.9T SP <T < T SP . Hence 

J 2 

V ~ 0.2534 



T 2 
J SP 



q ~ 0.0228^- (36) 

1 SP 

As a check of the correctness of the results of the Ginzburg Landau expansion, we can compare the prediction for the 
specific heat jump from Eq. I|34ll . AC V — — 1-4(1)'7 w ^ n ^ ne va hie given by the TBA mean field theory in Eq. 
Q32|l AC V ~ 1.39^. The 1% agreement between these two values provides a confirmation of the correctness of our 
Ginzburg-Landau expansion. From the behavior of the gap, we can obtain the behavior of the magnetic correlation 
length timag(T). If the gap is A, the magnetic correlation length at T — is u/A. Thus, neglecting thermal effects, 
we would obtain £ maff (T) ~ J/Tgp(l — T/Tsf)^ 1 ^ 3 ■ Near the transition, this correlation length becomes much larger 
than the thermal correlation length t; t h — u l i^T). This means that the exponential decay of the magnetic correlation 
that would result from the gap A(T) is completely masked by the thermal fluctuations which lead to a much shorter 
correlation length. The above results are valid for the case of a uniform dimerization S. In reality, near the transition, 
the dimerization can vary with the spatial location. To take into account the energy cost of these fluctuations, a 
(V(5) 2 term must be included in the Landau-Ginzburg effective theory. The bosonization approach allows to calculate 
the coefficient of this gradient term as outlined in App. [S] The full Landau Ginzburg free energy is now given by 



F L = Jdx [|(V5) 2 + |(T - T SP )5 2 {x) + \5\x)] (37) 
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where the constant cq measures the rigidity of the order parameter 5 and is given (see App. [SJ by: 

9 /tt\1/2 f3(2) J 4 a 



Co 



ir 2 \2J T(3/4) 4 T 3 1 J 



where /3(2) ~ 0.91596 ... is Catalan's constant^ and p, q are given by l|36[l . It is interesting to note that the coefficient 
of the gradient term falls as T 3 . We note that the Ginzburg-Landau coefficients calculated in Ref . l68l using a fermionic 
renormalization group treatment have the same dependence on T$p as in the present mean-field calculation. However, 
we do not expect an agreement of the numerical prefactors as the model studied in Rcf. 68 is different from the one 
studied here. The structural correlation length close to the transition can be evaluated from H37[) and is found be to 
be 

e{T) = * =(V£-Vfl2) ; Tsp 



p{T - T S p) V ZkTsp ) T - T S p 

Near the transition, £ 3> £th- This justifies the Landau-Ginzburg approach where the magnetic fluctuations are 
integrated out and only structural fluctuations close to the transition are retained. As was done in Re f. I68L the 
contributions of these structural fluctuations to the specific heat can be analyzed by the techniques of Refs[ 



D. Magnetic susceptibility 

Here we consider the effect of a magnetic field on the spin-Peierls system. The field can close the spin triplet gap 
and induce incommensuration. Here, we restrict ourselves to fields much smaller than the gap and study their effect 
on the spin Peierls transition temperature and the susceptibility. Using the perturbative approach^ generalizing the 
one of Sec. II Bl we find that for small magnetic fields there is a reduction of the transition temperature Tsp i.e., 
Tsp{h) = Tgp(0) — \h 2 , with A > 0. Using bosonization (see App.[SJ, we find: 

A-ffiU "* (40) 
ir z ±sp lbir^lsp 

This result is in accord with that of Ref. Ffol where it was found that A ~ 14.4/(167r 2 Tsp). For large fields, a similar 
calculation can be done provided that the field is smaller than the soliton gap and the system does not exhibit a 
transition to the incommensurate phase. 

On the other hand to calculate the finite temperature susceptibility, we need to generalize the TBA equations to 
include the effect of a magnetic field. To recapitulate, in the absence of a field, wc have two solitons of spin ±1, a 
light breather of spin forming a triplet and a heavy breather of spin 0. A field breaks the spin degeneracy and 
for small enough fields which do not induce any incommensuration wc can use the TBA to calculate the magnetic 
susceptibility. Following Ref. 0, we obtain the following TBA equations: 



eJff) = -M cosh (9- ha + — V 

2 9.TT ^ 



T 
2^ 



<T = -1,0,1 



de'K 11 (e-e , )ii-i(i + e- e " ( - e '^ T ) + — I de'K 12 {e-e')hx{i + e-^ 1 "'' T \ 



U0) 



Mcosh9+— V ( d6'K 12 {9-6')\n{l + e- l ^ 'V f ) + — [ dd'K 22 (9 - 6') ln(l + e "«( e ')/T ) 

„ — i n i ** « 



~Y 

K T 
2J 2tt 



dOM cosh 6 



J2 hi(l + e-^W/ ? ) + \/31n(l + e- 

.<T=1,0,-1 



M 
8%/3 



(41) 



where the € a {0) denote the reduced pseudoenergies of the solitons (<J = 1), antisolitons (a = —1) and light breathers 
(<t = 0) and h = h/J. It is easy to see from the equations above that e„(6) — eo(0) — her. This allows us to reduce 
the set of TBA equations to two: 
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^Mcosh0 + f^ / dO'Kuie - 0')ln(l + e"^')^) + ±- I d0'K 12 (6 - 9') ln(l + e"^'^) 

Z7T 



T 
2^ 



sinh" 



( 2 t) 



cosh' 



2 Z7T 



-s 2 (e')/f 



f d6'K 12 {6- 6') In 
2tt J 



1 



smh 2 ( A 



cosh 2 (^) 

In the presence of a magnetic field, the law of corresponding states l|28|) now reads: 



(42) 



T 2 



J 



T 



Tsp Tsp 



The magnetic susceptibility, x(T) = — lim^o^r?, satisfies the scaling relation 



1 



X(T) = -J% 



( T 



J v \T S i 



,0 , 



(43) 



(44) 



where T' y ' — dyJ r (x,y). For T > T$p, the susceptibility is that of a free Bose gas: x{T) = -^tj- For temperatures 
< T < Tsp, the susceptibility can be obtained numerically from l)42|l and l|44[l. The results are plotted in Fig. 0] 

As in the case of zero magnetic field, one has an effective field dependent Landau Ginzburg functional which 
describes the physics in the vicinity of the transition. For T < T$p the behavior of the magnetic susceptibility is 
obtained from the following Landau- Ginzburg expansion: 



F(T,h,8) = $(T-Tsp{h))6* + ±S 



q xi X0.2 



(45) 



2 V " v " 4 2 ' 
where xo — ^tj is the susceptibility of the undistorted chain. Minimizing F with respect to S, one finds for T < T$p 



The definition of \ then gives 



F(T, h) = -|-(T - T SP (0) + A/i 2 ) 2 - ^-h 2 . 



X (T < T SP ) = Xo + — (T - T SP (0)). 
9 



(46) 



(47) 



This behavior of the susceptibility is reminiscent of the one seen in a Neel antiferromagneliS which stems from the 
similarity of the mean field equations for the spin Peierls problem and the quasi-lD antiferromagnet. Using Eq. 14U|) 
and Eqs. i|36[) it is easily seen that the resulting susceptibility H47|) satisfies to the law of corresponding states 144fl . 

Numerically, one finds that = j'^ 7 p by fitting the susceptibility calculated near the transition to the Landau- 
Ginzburg form. If on the other hand, we use the values of p and q obtained in Sec. Ill Cl combined with the value of 
A in Eq. (|4*0"|) we see that value ^ = ^%ijt'J p ~ jTsp differs from that obtained in the presence of a field by less 
than a percent. A Keesom-Ehrenfest relation^ exists between the jump of the specific heat and the jump in the slope 
of the magnetic susceptibility as a function of temperature. This Keesom-Ehrenfest relation is given by the Landau 
Ginzburg theory as: 



rp2 (IT 

Tsp Aa, 



% tt 2 -^)(3/4) 2/3(2) 



4tt 2 



0.3724 



(48) 



Such a proportionality has been observed in experiments on TTFAuBDT in magnetic fields^ 
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FIG. 4: The magnetic susceptibility x(T) versus the reduced temperature; for T > Tsp, x(T) = l/(ir 2 J). 



E. Low temperature expansions 



In the preceding sections, we have obtained analytical results for T = and for T < T$p- In fact, the TBA 
equations are amenable to analytical study for T > (more precisely < T <C A(T = 0)). In this regime, we expect 
that the mean field gap A(T) remains very close to A(T = 0), so that the thermodynamics does not differ from the 
one of the Heisenberg chain with dimerization. With this assumption, the TBA equations (for ft, = 0) to lowest order 
are: 



d(0) = A(T = O)cosh0 + O(e 
e 2 (0) = A(T = O)V3cosh0 + O( 



A(T = 0) , 



(49) 



Substituting these in 11221) for the free energy F, we see that the correction to F is indeed O(e~ A ' T=0 - ) / T ) which 
justifies our original assumption that A(T) ~ A(0). We now derive low temperature expansions of the various 
physical quantities. It is convenient to use the the zero temperature dimerization <5o — 8(T — 0), the zero temperature 
gap Ao = A(T = 0), and the zero temperature groundstate energy Eq, the expressions of which are given in J^J, to 
express the corresponding finite temperature quantities as a function of the dimensionless variable 8. We obviously 
have S(T) = 6 Q S(T), and A(T) = A 5 2 / 3 (T). The total energy at T = reads 



E(S) 



4 /r(i/6)V 

V3\T(2/3)J 



T(3/4) 



r(i/4) 



7T 2 V2/ 



N 1/4" 


4 J 3 


72 ^ ? 4 

^ (5 3 




k 2 


2 



and it is easy to see that E(5) has a minimum for 8 — 1. Expanding around this minimum we find: 

2 



E{8) - E a 



A 2 



(s-iy 



(50) 



(51) 



The expression (|51|l is not the full expression of the free energy for T > 0, as we have also to take into account the 
contributions of the solitons and the breathers that are thermally excited. Since the heavy breathers have mass 
as can be seen from l|49|) . their contribution at low temperature is negligible with respect to the soliton contribution. 
Therefore, to lowest order, the thermal contribution to the free energy reads: 



Fsol. = 



Ao 
T 



^2/3 



(52) 



where K\ is a modified Bessel function*!* so that the full variational free energy is: 



(53) 
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Minimizing (|53|l with respect to <5 we obtain: 

7 , 9V3 



To lowest order, this gives: 



and: 



T 



8=1- 3 V < 



T5-V 3 
A 



2T 



A(T) = A 



ttA 



(54) 



(55) 



(56) 



Substituting l|55(l in (|53(l . we see that the correction to the elastic energy plus ground state energy of the dimerized 
chain is of order e - 2A o/ T is is therefore negligible compared to the contribution of the solitons. In physical terms, 
this means that at sufficiently low temperature, the thermodynamics of the spin-Peierls chain is the same as the 
thermodynamics of a chain with a constant dimerization. Using this result, we find that: 



A 2 

F(T) = 



^Ao Kl It" 



which leads to a low temperature specific heat of the form: 



C V (T) 



3yg A ( A 

71- 3 / 2 J 



T 



3/2 



e T 



o(T 



-3/2„-A /T 



(57) 



(58) 



In the presence of an infinitesimal applied magnetic field h <C T, the lowest order contribution to the low temperature 
free energy is: 



F(T) = - 



A 2 



IT 



Ao 
T 



(1 + 2 cosh-) 



(59) 



and the magnetic susceptibility is readily obtained as: 



7T Z J 



(60) 



III. CONCLUSIONS 



In this paper, we have studied the thermodynamics of the spin-Peierls system treated within a mean-field approx- 
imation. Using a combination of bosonization methods and the thermodynamic Bethe ansatz, we have been able to 
obtain quantitative results for the spin-Peierls transition temperature T$p, the spin-Peierls gap to triplet excitations, 
the specific heat and magnetic susceptibility at arbitrary temperatures. Our calculations are a quantitative improve- 
ment of the results obtained by Cross and Fisher (who were restricted to the vicinity of the spin-Peierls transition 
temperature) and consequently help us obtain the effective Landau Ginzburg functional that describes the physics of 
dimerization close to the transition. It would be interesting to study this Landau Ginzburg theory in one dimension, 
following Ref. to understand more quantitatively how lattice fluctuations affect the thermodynamics. Similarly 
to Ref 68, we should expect a regime of renormalized Gaussian fluctuations for 0ATgp F < T < TgjF ', and a regime 
dominated by kinks for 0.3Tsp < T < OATsp as long as one dimensional fluctuations dominate. Also, it should be 
possible to use the Landau- Ginzburg description to study the three-dimensional ordering of the dimerization along 
the lines of Refs. 67 69. The dependence of the transition temperature on interchain coupling will be similar to the 
one predicted in Ref. l74l since both models belong to the same universality class. A more direct extension of the 
present work would be to study the commensurate-incommensurate transition driven by an external magnetic field 
and then comparing the predicted results to various experiments on spin-Peierls systems. It would also be interesting 
to extend this study to the generalized spin-Peierls transition obtained in ladders under magnetic field 60 or to the 
antiferromagnetic phase transition obtained in the same system^ These questions are left for future work. 
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APPENDIX A: CALCULATION OF THE RIGIDITY 

Here, we present a derivation of the rigidity in our Landau-Ginzburg effective action. In the continuum limit, the 
space dependent dimerization leads to the modification of the sine Gordon term in J2J 



Hi, 



2.9 



(2-Ka) 1 



dx5{x) cos y/2<j>. 



(Al) 



Close to the transition, the second order correction to the free energy of the spin chain induced by the spin-phonon 
coupling is given by: 



F s 



1 7ra / g \ 2 



4 (3u \Tra 



(-)' 

\7ra/ 



dx' I S(x)S(x')x(x — x' , r), 

-oo JO 



and: 



x(x - x', t) = \f2 



2nx 2ttt 
cosh — cos — — 

pu [3 



(A2) 



(A3) 



In Fourier space, Fg oc J ^8(q)d(—q)x(q,i^ n = 0). To obtain the gradient term, it thus suffices to calculate the 
Fourier transform \- In the limit q — > 0, ic{q,uo n = 0) = x(0, 0) + g 2 /2x"(0, 0). Plugging this form into l|A2(l . it is 
straightforward to find the rigidity. 

To find x(q) we generalize slightly the calculation of the q = response function, and consider: 



X(q) = J dxdT X (x,T)e lqx 
= J dxdr 



V2 



cosh 

Using Eq. (8.12.5) in Ref. 153. we can rewrite the integral: 

dv 



(A4) 



cosh v — cos cosh 



1 / 2ttt 



(A5) 



Consequently, to calculate (|A4|) . we only need the integral: 



P p (_ 2ttt\ _ n£ 

'o -w^l co& P r r -[r(i + ^)r(l-i^)p : 

which is easily obtained from Eq. (8.14.16) in Ref. l52l The final result is: 

x(q) = 9 9 ■ 

2cosh(^)r(| + ^) 2 r(|-i^) 2 

Expanding x{l) to second order in q, we find: 



F = F n 



dq 



.9(9)| 2 



where 



167r 2 a 3 Tr(3/4) 4 J 2tt 

g(q) = [ gS(x)e^ x dx, 



(A6) 



(A7) 



(A8) 



(A9) 
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and 7] = * ~%J 3/4} and ^(x) is the tri gamma function (see Ref.H3p.260). The number V> (1) (3/4) can be expressed 
as a function of Catalan's constant 5 — (3(2) as ip^ft/A) = tt 2 — 8/3(2) ~ 2.5419. Finally, using the expression of g(x) 
as a function of S(x) we obtain the rigidity cq as: 



Co 



64 



— (-) 

8tt 2 \2/ 



^W(3/4) 



1/2 j3(2) 



TT 

2 
J 4 a 



1/2 



1 J 4 a 
T(3/4) 4 Y T 



r(3/4) 4 r 3 



(A10) 



A similar calculation can be done to obtain the reduction of the critical temperature as a function of the magnetic 
field. When the system is magnetized, incommensuration in the staggered operator sets in as the wavevector shifts 
from ir/2 to ir/2±h/u and the equation giving the critical temperature reads: 



2 \2anJ 



/3sech(^) 



r r (! + <£) 



(All) 



This implies that: 



T SP (h) cosh 



2T SP (h) 



iwT SP (h) 



4 4TrT SP (h) 



T SP (h = 0)r(3/4) 4 



(A12) 



the Equation l|A12(> was obtained in Ref. [7CJ using a real time calculation of the response function. This can be seen 
explicitly by expressing the infinite products in Ref. IJiJin terms of Gamma functions. Expanding IJA12J) around small 
h, one obtains for magnetic fields h <§; Tsp(0): 



T S p(h) 



1-2(tt 2 -V> (1) (3/4)) 



T S p(0) _ v ' T y "' '" \4ttTsp 

In terms of Catalan's constant, the spin-Peierls transition temperature in the presence of a small field is given by 



(A13) 



T S p(h) 
Tsp(0) 



1-/3(2) 



(— 

\ttTsp 



o(h 2 ). 



(A14) 
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